Decoherence and thermalization of a pure quantum state 
in quantum field theory 
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We study the real-time evolution of a self-interacting O(N) scalar field initially prepared in a pure, 
coherent quantum state. We present a complete solution of the nonequilibrium quantum dynamics 
from a 1/JV-expansion of the two-particle-irreducible effective action at next-to- leading order, which 
includes scattering and memory effects. We demonstrate that, restricting one's attention (or ability 
to measure) to a subset of the infinite hierarchy of correlation functions, one observes an effective 
loss of putity /coherence and, on longer time scales, thermalization. We point out that the physics 
of decoherence is well described by classical statistical field theory. 
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Quantum decoherence is a fundamental process whose 
understanding is a central issue in many areas of physics. 
Topical examples include measurement theory in quan- 
tum mechanics jll], the physics of Bosc-Einstein conden- 
sates 0, or of quantum computers 0, neutrino physics 
4 1, high-energy nuclear collisions Q, black hole physics 
6], or the description of primordial fluctuations in infla- 
tionary cosmology @, Q ■ It is a genuine nonequilibrium 
process, which requires the real-time description of quan- 
tum dynamics. Analytic descriptions can be obtained for 
exactly solvable models and/or simple enough approxi- 
mations, e.g. assuming a linear coupling between the sys- 
tem and its environment, neglecting back-reaction, etc. 
A complete microscopic description in realistic quantum 
field theories (QFT) is a notoriously difficult task |9j, 
which requires first principle calculations of the nonequi- 
librium quantum dynamics. 

Recent years have witnessed substantial progress con- 
cerning the description of quantum fields out of equilib- 
rium [lOj ■ Two- particle-irreducible (2PI) functional tech- 
niques provide a powerful tool to devise systematic, prac- 
ticable approximation schemes, valid for arbitrarily far 
from equilibrium situations [H, flo| . It has been demon- 
strated that a coupling or 1 /TV-expansion of the 2PI ef- 
fective action at lowest nontrivial order can describe far- 
from-equilibrium dynamics and subsequent (quantum) 
thermalization without further assumption |lll4l3| |. In 
this letter, we use 2PI techniques to compute the dy- 
namics of decoherence from first principles in QFT [22j . 

There are various uses of the concept of (de) coherence 
16], the most widely discussed being the so-called 
environment-induced decoherence, which results from the 
interaction of the system under study with a (thermal) 
bath of unobserved degrees of freedom (d.o.f.) [l[. Even 
in the absence of an environment, decoherence of a sub- 



*APC is unite mixte de recherche UMR7164 (CNRS, Universite 
Paris 7, CEA, Observatoire de Paris) 



set of d.o.f. may result from some kind of coarse-graining 
[l7l [HI . Here, we adopt a different point of view, also 
advocated e.g. in || 0, [TH, Il9j . Even if one keeps all 
the dynamical d.o.f., reconstructing the actual state of a 
given (closed) system requires a precise knowledge of its 
independent correlation functions. In practice, however, 
such information is often not experimentally accessible 
and one has to infer the state of the system from the sub- 
set of measured correlation functions. This "incomplete 
description" picture actually underlies the very concept 
of thermalization [| 01]. Similarly, a system prepared in 
a pure quantum state may appear as a statistical mixture 
to the observer who has only partial information. 

In this letter, we show that, starting from a pure, co- 
herent quantum (Gaussian) state, an observer who only 
measures the subset of equal-time two-point functions, 
observes an effective loss of purity/coherence and, even- 
tually, (apparent) thermalization. We study a relativistic 
self-interacting O(N) scalar field and present a complete 
numerical solution of the nonequilibrium dynamics from 
a 1/iV-expansion of the 2PI effective action at next-to- 
leading order (NLO) [l3j]. The approach allows us to 
study the strong coupling regime and is a valid descrip- 
tion for states with a high degree of quantum coherence 
which, as explained below, are characterized by strong - 
possibly non-perturbative - field fluctuations. 

A pure quantum state remains such under a uni- 
tary evolution. Here, although the complete information 
about the initial Gaussian state is contained in the set of 
two-point functions, nontrivial higher correlators develop 
in time due to the non-Gaussian dynamics: The informa- 
tion about the initial state spreads in the space of corre- 
lation functions and gets lost to our observer, resulting in 
the effective loss of purity/coherence and, at late times, 
the thermalization of the effective density matrix. We 
stress that no environment of incoherent d.o.f. is needed. 
We illustrate this point by preparing a completely pure 
initial state. Finally, we show that the regime of de- 
coherence, characterized by strong field fluctuations, is 
well-described by classical (statistical) field theory. 



We consider a relativistic real scalar field <p a (a = 
1, . . . , N) with classical action 

<%] = -J d 4 x |i(^ a (D + m 2 )( y 5 a + ^ : (<^ a ) 2 j (1) 

where summation over repeated indices is implied. For 
vanishing field expectation value, (tp a (x)) = 0, the corre- 
lation functions of the quantum theory can be obtained 
from the 2PI effective action T[G], parametrized by 
the time-ordered connected propagator (Tip a (x)tpb(y)) = 
S ab G(x,y) 0. The 1 /iV-expansion of T[G] at NLO in- 
cludes the infinite series of diagrams shown in Fig. [TJ 
where the lines represent the propagator G fl3| . 

To set up the initial-value problem, it is useful to de- 
compose the two-point function G into a statistical (F) 
and a spectral (p) component, both real functions: 

G(x, y) = F(x, y) - ^sign c (x° -y°)p(x,y). (2) 

Notice that F(x,y) = F(y,x) and p{x,y) = —p(y,x). 
For Gaussian initial conditions, the equations of motion, 
obtained as 5T[G]/SG(x, y) = 0, read 

r° 

[D x + M 2 (x)]F{x,y) = - / d i zZ p (x,z)F(z,y) 

Jo 



p x + M 2 (x)] p{x,y) 



d i zY, F (x,z)p(z,y),(3) 



d 4 zT, p {x,z)p(z,y). (4) 



where J Q X d i z = J* dz° Jdz. The effective mass term is 
given by 



^ ^ ^ N + 2 F(x,x) 



M 2 (x) = m 2 + A- 



67V 



(5) 



and the self-energies by [l3( 
^_ 

A 

3N 



E F (x,y) 
^ P {x,y) 



F(x, y)I F (x, y) - -p(x, y)I p {x, y)J , (6) 
p(x,y)I F (x,y) + F(x,y)I p (x,y) . (7) 



The functions I F and I p resum the infinite series of bub- 
ble diagrams in Fig. [T] 



If(x,v) = U F (x,y) 



l P (x,y) = K P (x,y) 

with the elementary bubble 
A 



d 4 zl p (x, z)U F (z,y) 



d 4 zI F (x,z)U p (z,y), 



d 4 zl p (x, z)U p (z, y), 



Tl F (x,y) 



F (x,y) - |P (x,y) 

F(x,y)p(x,y). 



(8) 
(9) 

(10) 
(11) 




FIG. 1: The dots indicate that each diagram of the series is 
obtained from the previous one by adding another bubble. 

The self-energy kernels E F and E p describe scattering 
and memory effects and are responsible for non-Gaussian 
correlations to develop in time, a key ingredient in the 
present discussion. Gaussian - collisionless - evolution 
equations, such as leading-order Large-iV, or Hartree 
approximations, where the RHSs of Eqs. (©-(j!]) vanish 
identically, are known to fail to describe phenomena such 
as damping of unequal-time correlators or thermalization 
[lOl ]. We show below that they also fail to capture the 
physics of decoherence in the present approach. 

We consider spatially homogeneous and isotropic 
states, for which F(x,y) = F(x°,y°, |x — y|). Accord- 
ingly, we introduce the Fourier decomposition (p = |p|) 



F(t,t', 



dp 

(2tt) 3 



■F p (t,t'), 



(12) 



and similarly for p(x,y). We consider an observer who 
only has access to the subset of equal-time two-point 
functions: F p [t) = (^(i)ft>(*)>» #p(*) = + 
TT^(t)if p (t)) and K p (t) — (7r£ (t)7T p (t)). The mimimum- 
bias state compatible with these measured observablcs 
is described by an effective Gaussian density matrix [l9| 
D e fi(t) — rip-^pWi where the product runs over inde- 
pendent d.o.f. in momentum space with, up to a normal- 
ization, 

D p (t) ex exp {n p {t) [F p (t)nl (t)ir v (t) + K p (t)^ p (t)^ p (t) 
-R p (t)(7rl(t)<p p (t) + 4(t)7r p (t))] } , (13) 

where K p (t) = — ln[l + 1 / n p (t)] / (2n p (t) + 1), with 

n P (t) + \ = y/F p (t)K p (t) - B? p {t) = ap(t) . (14) 

Contrarily to the correlators F p (t), R p (t) and K p (t), 
which can be modified by a canonical redefinition of the 
field variables (if, it), n p (t) is a canonical invariant and, 
actually, the only truly intrinsic property of the density 
matrix (fT3|) [8]. It provides an absolute measure of the 
quantum purity of the system's state through 23] 



tr[D»(t)] = (2n p (t) 



1 



< 1 



(15) 



which equals 1 for a pure state. Note also that, when- 
ever the system is well-described by a thermal ensemble 
of weakly interacting quasi-particle (QP), n p (t) defines 
a QP occupation number and follows a Bose-Einstein 
distribution as a function of the QP energy 13]. We 
shall make use of this property to characterize the de- 
gree of thermalization of the system, although we stress 
that none of our results rely on such a QP interpretation. 
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FIG. 2: The function ln[l + l/n p (t)] as a function of e p (t) for 
times Mot = 2.5x2" for n = 0, . . . , 14. A straight line at late 
times corresponds to a Bose-Einstein distribution. The inset 
shows the rapid decrease of tr[£)p(t)] as a function of time for 
modes (from bottom to top) p/a = 1, 5, 10, 15, 20, 24 where 
a/Mo = 0.15. Higher momentum modes, p/a = 25, 27, 30, 
are also shown for comparison. 

Note finally that n„(t) is time-independent for Gaus- 
sian approximations |24j . Therefore, no loss of quantum 
purity, in the sense discussed here, occurs for such ap- 
proximations. We stress that this is not in conflict with 
existing studies of decoherence in free field [7| or mean- 
field approximations [l8[, which typically consider an ef- 
fective coarse-grained Gaussian density matrix where one 
averages out some rapidly varying d.o.f. This results in 
an effective loss of information and hence of quantum 
purity /coherence, due e.g. to dephasing [l8[, already at 
the Gaussian level. In the present incomplete descrip- 
tion approach, decoherence is due to the non-Gaussian 
dynamics, i.e. direct scattering and memory effects [2o| . 

We now come to the discussion of quantum coherence 
which, unlike purity, is a basis-dependent notion. The 
equal-time two-point correlators can be parametrized as 



e p {t)F p {t) 

Rp(t) 

K p (t)/e p (t) 



a p (t)[l - 7 P (t) cos0 p (f)] 
-a p (t)j p (t) sm(f> p (t) , 
a p (t)[l +l P (t) cos</> p (£)] 



where 



a p (t) 



K p (t)+el(t)F p (t) 
2e p {t) 



n p (t) 



(16) 
(17) 
(18) 



(19) 



provides an alternative definition of a QP occupation 
number n p (t). Note that < n p {t) < n p (t). Here, the 
energy scale e p (t) defines the QP basis in which to dis- 
cuss decoherence. We use e p (t) = \/p 2 + M 2 (t), which 
we found gives a good description of the oscillation fre- 
quency of two-point correlators. The basis-dependent oc- 
cupation number n p (t) is related to the canonical invari- 
ant n p (t) through the coherence parameter [26| 



7p (t) = yjl-al(t)/al(t). 



(20) 



The case 7 p (i) = or, equivalently, n p (t) = n p (t), cor- 
responds to the thermal-like density matrix D p (t) oc 
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FIG. 3: Exponential decay of the coherence parameter 7 P (t) 
as a function of rescaled time r = (A e ///6) 2 (Moi) (see text), 
for the mode p/Mo = 0.15, for four couples of parameters 
(70, A): (0.99,3.05) [solid] and (0.95,6.62) [long-dashed], for 
which A e// = 21.7, (0.95,2.2) [dashed], for which \ eff = 
7.215 and (0.95,0.915) [dotted], for which X eff = 3.0. The 
inset shows the corresponding decay rates, obtained from ex- 
ponential fits, for modes p < p c . The agreement of the first 
two runs illustrates the A e //-scaling of early-time decoher- 
ence. The overall agreement shows the approximate \&ff- 
dependence of the decoherence rate. 

exp{« p (i)F p (i) [nl(t)ir p (t) + e 2 p (t)<pl(t)v P (t)]}. This in- 
cludes the vacuum-like state: n p (t) = -f p (t) = 0. More- 
over, 7 P (i) controls the size of off-diagonal matrix ele- 
ments of the density matrix in the so-called two-mode co- 
herent state basis [27| : The latter exhibits non-trivial cor- 
relations (i.e. quantum coherence) between macroscop- 
ically distant semi-classical states for j p (t) —> 1. Note 
that this limit is also characterized by strong field fluc- 
tuations since the correlators F p (t), R p (t) and K p (t) are 

ocl/Vl-7p 2 (*)>!• 

At time t = 0, we prepare modes with p < p c in a 
pure, highly coherent quantum state characterized by 
n p (0) < 1, 7 P (0) = 70-1 and cf> p (0) = tt/2, while 
modes with p c < p < A, with A the ultra-violet cut- 
off, are prepared in a vacuum-like state: n p (0) ^ 1 and 
7 P (0) = [28|. We emphasize that there are no incoher- 
ent d.o.f.: The initial state is a completely pure quan- 
tum state. Eqs. (|3^- ([Tl~T) are solved numerically with- 
out further approximation (except for discretization of 
time/momentum- integrals and truncation of memory- in- 
tegrals). We show results for p c /M — 3.6 and A/M = 
10, where Mo = M(0) is the initial effective mass. 

Fig. [2] presents an overview of the time-evolution. It 
shows snapshots of the function ln[l + l/n p (t)] as a func- 
tion of e p (t) at various times. One observes a substan- 
tial growth of n p (t) for all modes, signaling the effective 
loss of quantum purity at early times, followed by a slow 
approach to an effective quantum thermal equilibrium, 
characterized by a Bose-Einstein distribution. The loss of 
quantum purity is further illustrated in the inset, which 
shows the rapid decay of tr[Dp(i)] for modes p < p c . 

After the loss of quantum purity, we find that n p {t) ~ 
n p (t), signaling a corresponding loss of quantum coher- 
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ence according to Eq. (|20|) . Fig. [3] shows the time evo- 
lution of the coherence parameter 7 p (t) for the mode 
p/Mo = 0.15, for various sets of parameters. In all cases 
we find that decoherence is well-described by an exponen- 
tial law. The inset shows the corresponding decoherence 
rates for the modes p < p c - 

Classical scaling: In the highly coherent limit, 70 — > 1, 
one has, roughly speaking F ~ 1 / \J 1 — 70 ^> p ~ 1 
which signals the enhancement of classical vs. quan- 
tum fluctuations [l2| . In this regime one can neglect the 
second term in brackets on the LHSs of Eqs. (fTU|) and 
©. It is then easy to check that under the simultaneous 
rescaling of the correlators and the coupling constant: 
F — > r)F, p — s> p and A — » A/77, with rj an arbitrary con- 
stant, the F and p-components of II, / and £ scale just as 
F and p respectively and that Eqs. ©-(HI) are left invari- 
ant. This is characteristic of the regime of strong field 
fluctuations, where a classical (statistical) field theory de- 
scription is appropriate. Indeed, in classical field theory, 
the action being defined up to a multiplicative constant, 
a rescaling of the field (i.e. of the initial conditions) can 
be entirely absorbed in a change of coupling. Therefore, 
as long as the system is highly coherent (^f p (t) ~ 1), we 
expect the dynamics not to depend separately on 70 and 
A, but instead on the combination A e // = A/77 where 
V = yl — 7o ■ We checked, from our numerical simula- 
tion of the full quantum dynamics, that this is indeed the 
case during the early-time decoherence regime, as illus- 



trated in Fig. [3) Runs with different values of 70 and A 
but the same A e // are essentially indistinguishable. Fur- 
thermore, we observe that, despite the rather strong (ef- 
fective) couplings employed in some simulations, the de- 
coherence rates approximately follow a perturbative-like 
Agyj-scaling, as also illustrated on Fig. [3] 

As a final remark, we mention that fixing the value 
of the initial mass Mo absorbs a large A 2 -dependence. 
This simple, though approximate, renormalization turns 
out to be numerically sufficient for the results presented 
here. We find that, if the late time thermalization is 
cutoff-dependent, a known artifact of Gaussian initial 
conditions [2l|, the regime of effective loss of quantum 
purity/coherence is largely cut-off insensitive. 

To the best of our knowledge, this work provides the 
first complete microscopic description of the process of 
decoherence in a realistic QFT. It is an exciting obser- 
vation that the relevant dynamics is well-described by 
classical statistical field theory, which can be solved ex- 
actly by means of standard Monte-Carlo Techniques. A 
particularly interesting question is to investigate how the 
present results are modified as one includes the knowl- 
edge of higher order equal-time correlation functions in 
the set of measured observables. 
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